Log-Poisson Cascade Description of Turbulent Velocity Gradient Statistics 
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The Log-Poisson phenomenological description of the turbulent energy cascade is evoked to discuss 
high-order statistics of velocity derivatives and the mapping between their probability distribution 
functions at different Reynolds numbers. The striking confirmation of theoretical predictions sug- 
gests that numerical solutions of the flow, obtained at low/moderate Reynolds numbers can play 
an important quantitative role in the analysis of experimental high Reynolds number phenomena, 
where small scales fluctuations are in general inaccessible from direct numerical simulations. 
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I. INTRODUCTION 

Since the pioneering experimental work of Batchelor 
and Townsend, published exactly sixty years ago [Tl| it 
is known that scale dependent galilean invariant observ- 
ables, like velocity differences, fluctuate in a strongly non- 
gaussian way at small scales in turbulent flows. This kind 
of statistical behavior, generally referred to as "inter- 
mittency" , indicates that the K41 picture of turbulence 
0, Q , which actually would correspond to the existence 
of a homogeneously distributed energy dissipation field 
should break down, a fact notoriously anticipated 
by Landau as early as in 1942 Q. Not less remark- 
ably, long before additional breakthrough experiments 
were performed phenomenological models of the en- 
ergy cascade advanced the conjecture that intermittency 
should be related to the stochastic multiplicative nature 
of the energy cascade process 0, Q , implying that small 
scale strong fluctuations are, in some sense, fed by the 
weaker large scale ones. 

The intermittency phenomenon is commonly associ- 
ated with the anomalous scaling of velocity structure 
functions. A comprehensive description dealing with 
both anomalous scaling and the non-gaussian behavior 
of intermittent observables is a major challenge of three- 
dimensional turbulence theory [!||. Small scale strong 
fluctuations are believed to reflect the dynamics of coher- 
ent structures like vortex filaments. Even though this is 
a very open problem, a similar physical picture is in fact 
well-established in simpler contexts, as in Burgers turbu- 
lence , with shocks playing the role of "vortices" . 

The log-Poisson model [I4I yields perhaps the most 
intriguing description of the turbulent multiplicative cas- 
cade, since, as it is well-known, it leads to the accurate 
She-Leveque intermittency exponents of velocity struc- 
ture functions [13]. The phenomenological work of She 
and Leveque is also of great physical appeal, once it 
places vortex filaments as a fundamental ingredient in 
the production of intermittency. 

We are interested, in this work, to know what the log- 



Poisson model may tell us about the profiles of velocity 
gradient pdfs. We deal here with two sets of pdfs for fiows 
associated to different Reynolds numbers. One of them 
is obtained from an atmospheric surface layer experiment 
[13, 111] and the other from a direct numerical sim- 
ulation (DNS) of homogeneous and isotropic turbulence 
[13]. The underlying motivation in this choice of sys- 
tems is to show that numerical low/moderate Reynolds 
number results can be useful in the modelling of fiows 
that cannot be directly simulated (even in a foreseeable 
future). The very same claim was put forward in a pre- 
vious letter [ll], where, despite the force of evidence, 
lacked some phenomenological basis, which, then, we de- 
velop here. We find that a bridge between low and high 
Reynolds number pdfs can be built within the framework 
of the log-Poisson model [l^ . 

This paper is organized as follows. In Sec. II we briefly 
review the multiplicative cascade models, introduce the 
log-Poisson model and compute hyperflatness factors of 
velocity gradient fluctuations, comparing them to recent 
estimates. Two relevant theorems related to velocity gra- 
dient pdfs are also established. In Sec. Ill, we present 
the experimental and numerical data that was analysed. 
In Sec. IV, the experimental and the numerical velocity 
gradients are closely matched with the help of a Monte- 
Carlo procedure based on the theorems of Sec. II. In Sec. 
V, we summarize our results and point out directions of 
further research. 



II. VELOCITY GRADIENT STATISTICS 

Multiplicative Cascade Models 

In the multiplicative cascade models [3], one assumes 
that energy flows from the integral scale L to the dissipa- 
tive scale 77 through a number of "quantum" steps asso- 
ciated to eddies of sizes L,_L/a, L/a^, where a > 1 is 
an arbitrary rescaling factor. At length scale Im = L/a™' 
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the fluctuating energy transfer rate is defined as 

Em = eoWlW2...Wm , 



(2.1) 



where the W^s are positive independent random vari- 
ables, with unit expectation value, (W) = 1, so that the 
mean energy transfer rate is conserved along the cascade 
process, i.e., (£,„) = eo. The scaling behavior of veloc- 
ity structure functions, Sq{r) = {(Sv)'^) is, then, derived 
with the help of Kolmogorov's refined similarity hypoth- 
esis, which postulates that fluctuations of 6v at scale 
have the same moments (up to constant numerical fac- 
tors) as (em^m)^/^. 

Analogous phenomenological arguments can be put 
forward to deal with the case of velocity derivatives - 
generically denoted in the following by dv. The essential 
idea is to assume that spatial fluctuations of the velocity 
field are smooth at the dissipative scale, and, therefore. 



(2.2) 



where, above, Sv,^ is the velocity increment defined at 
length scale ry. One may write, based on purely dimen- 
sional grounds, r] ~ {v'^ /eriY^'^ . Thus, substituting the 
latter on (|2.2p . we get 



dv 



(2.3) 



a statistical correspondence not unknown to the previ- 
ous literature 20] . A more interesting formulation of the 
refined similarity hypothesis is given in terms of proba- 
bility distributions. As it is clear, velocity gradient pdfs 
can be always written as 



p{dv) 



depi{e)p2{dv\e) , 



(2.4) 



where p2{dv\e) is the velocity gradient pdf conditioned 
on the energy transfer rate = e and pi(e) is the pdf 
associated to events which have = e. The refined 
similarity hypothesis is, then, the statement that at large 
Reynolds numbers. 



approximation, F{x) turns to have a gaussian profile, but 
one expects asymmetric corrections to be relevant in the 
problem of longitudinal structure functions, due to their 
non- vanishing skewness. 

Log-Poisson Model 

In the log-Poisson model [ll|, one writes down the 
energy transfer rate factors as 

W = a^~™ , (2.8) 

where a — 3/2, p = 2/3 and m > is a Poisson random 
variable, with expectation value 

'3' 



m a = 2 In , 

a- 1 V2 



(2.9) 



In order to cope with velocity gradient fluctuations, 
it is necessary to set up in first place the total number 
N of cascade steps associated to the turbulent flow un- 
der scrutiny. In other words, we would like to find TV, 
such that f] = L/a^ . We stress that the multiplicative 
cascade description addressed here is far from being a 
rigorous framework, since we take the Kolmogorov scale 
77 ^ to be a fluctuacting quantity. Thus, N should 

be defined, necessarily, from some averaging procedure. 
We adopt a simple prescription based on the definition 
of the Reynolds number as Q 



Re 



L3e = 



Therefore, we find 



N 



■ loga Re 



a3 



N 



(2.10) 



(2.11) 



An alternative and useful expression for N can be given 
in terms of the Taylor-based Reynolds number R\, which 
follows by taking the homogeneous isotropic result R\ = 



N 



6 



loga R^ 



■ log, 15 



(2.12) 



P2{dv\e) 



(2.5) 



where F{-) is a universal (Reynolds number independent) 
function of its argument. In fact, taking (|2.4p and (|2.5p . 
it is not difficult to show, in agreement with (|2.3[) . that 



{{dv)'^)^C,{{e,/:.)i 



where 



dxx^Fix) . 



(2.6) 



(2.7) 



It is worth noting that the form of the universal func- 
tions F{x) for the case of velocity differences has been 
the subject of experimental research [2ll,[23|. As a first 



Hyperflatness Factors 

As a direct application of the log-Poisson model, we 
compute the Reynolds-dependent velocity gradient hy- 
perflatness factors, defined as 



Hq{Rx) 



(2.13) 



A straightforward manipulation of (|2.13[) , taking into ac- 
count ([23]), jMl) and jlSl), gives 



H,iR,) = ^{('-^y)^^{Wir^A,Rl , 



(2.14) 
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where 



with 



j4g — |-15 



3 3 a/i . 



(2.15) 



(2.16) 



In particular, the skewness and flatness coefficients pre- 
dicted by (|2.16p are as ~ 0.13 and 0:4 = 1/3, respec- 
tively. Good support is found from the recent account of 
Ishihara et al. [23], which yields az = 0.11 ±0.01 and 
a4 = 0.34 ±0.03. 

If Ra and Rb are Taylor-based Reynolds numbers, re- 
spectively associated to flows with Na and Nb cascade 
steps, then (|2.14p implies that 



Hq(RA) 
HqiRs) 



and, thus, taking into account (|2.16p 

Na-Nb 



2a. 



3 , Hq{RA) 



Hq{RB) 



(2.17) 



(2.18) 



a quantity that measures the "distance" between cas- 
cades, going to play an important role in Sec. IV. 

Velocity Gradient PDFs 

We are interested to explore further consequences of 
the log-Poisson cascade picture in the setting of velocity 
gradient pdfs. In order to render the exposition more 
systematic, we introduce two important results in the 
form of theorems. 

Theorem 1. Let = {{dv)'^). The standardized pdf 
p{dv) = ap{adv) has a universal profile at fixed R\. 

Proof. We obtain, from ([^ and (^3]) . 

p{dv) — ap{adv) — 

00 

dtp\{t)\/ V I tF{a^ vjtdv) 

/ depi{v(j^e)^/TreF{./Tr€dv) . (2.19) 
^0 



Our task, thus, is to show that va"^ pi{v(T'^e) is indeed 
universal. Since the sum of Poisson random variables is 
also a Poisson random variable, Eqs. (12. ip and ()2.8p lead, 
for a cascade with N steps, to 



N fi—m 



(2.20) 



where m is a Poisson random variable with expectation 
value Nc. We may write, thus. 



^ (jVc)"'e 

m=0 



-Nc 



ml 



(2.21) 



Now, according to (|2.6p we write the variance of dv as 
— €1^02/ V, and, therefore, find 



ua^ pi{va^e) = C2 ^ 



(A^c)™e 



-Nc 



m=0 



which, in fact, ultimately depends only on R\. 



5{C2e - a^'^-'") , 

(2.22) 
□ 



Theorem 2. Let A and B denote flows with Taylor- 
based Reynolds numbers Ra and Rb, associated to log- 
Poisson cascades with Na and Nb steps, and velocity 
gradient pdfs 



PAidv) 



depti£)^/17reF{y^edv) 



PBidv) / depf{e)y^i'/eF{y/iy/edv) . (2.23) 



Lt follows that 

PA{dv) 



'-^Ki.)pBi-) , 

Q X X 

where K{x) is the pdf of the random variable 



X = 



^{NA~NB)ti- 



(2.24) 



(2.25) 



which, on its turn, is defined in terms of m, a random 
Poisson variable with expectation value {Na ~ Nb)c. 

Proof. A proof follows by direct substitution of the ex- 
plicit form of K{x) in l^^ . Defining g = a^''/^, with 
N = Na — Nb, we may write 



(iVc)'"e 



-Nc 



-6{x-ga ") 



(2.26) 



Using (p:^ and fTZE^ . we obtain, for the RHS of fTI^ . 

r '-^Kix)PBi-) 

Jq X x 



dx 



00 



E 

m=0 
.2 



(A^c)™e- 
m! 

.2 



-Nc 



-5{x-ga 2) 



X / devaj3p'({vcr%e)^/l/tF{y/l/edv/x) 
r dx ^ ^ (A^c)™e-^^ 







X 







m— ni'—O 



m 



{Nbc 



tin' -Nbc 



-5{x - ga~^)S{C2e - a^^^"") 

(2.27) 



m'! 

X ^/l/^F{y/T/~edv/x) . 

Performing the substitution e ^ e/x^ in (|2.27p and sub- 
sequently integrating over x, we get 

00 ^ 00^ (iVc)'"(iVBc)"'e-^^= 
^2 / de}^ }^ 



1 'm'' 



m— m'—O 



X S{C2e- g^a'^^''-"'-"'')y^eF{^/Tredv) 



(2.28) 
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Define, now, p = m + m', so that (|2.28p becomes 

(iVc)™(A^Bc)™-Pe-^^'= 



p—O m— 



m\[m — p)\ 
X 5{C2e- g'^a^^^-P)y/TreF{y/iredv) 

(2.29) 



p=0 

^i^(v/T7e9u) = pA{dv) 



□ 



In view of Tlieorem 2 we can devise a straiglitfor- 
ward Monte-Carlo integration procedure in order to re- 
late velocity-gradient pdfs defined at different Reynolds 
numbers. In fact, if x > and y are random variables 
of two independent stochastic process, described, respec- 
tively, by pdfs K{x) and p{y), then the random variable 
z = xy IS given by the pdf 



{5{z - xy))^^y 



dxdyK{x)p{y)5{z — xy) 

^ -K{x)p{-) = p{z) , (2.30) 



where we have used (|2.24p in the last equality above. We 
have found that it is greatly advantageous to use Monte- 
Carlo integration, instead of more traditional numerical 
methods, a fact probably due to the bad convergence 
properties of the latter in our particular problem. 



III. ATMOSPHERIC SURFACE LAYER 
EXPERIMENT 

Atmospheric surface layer velocity fluctuations were 
studied over a grass-covered flat surface in the Sils-Maria 
valley, Switzerland a place which hosts reasonably 
stable winds. The results reported in this work corre- 
spond to measurements of all of the nine components of 
the velocity gradient tensor, performed in a tower 3.0 
m high. The velocity signal was recorded at sampling 
rate of 10 KHz (which was high enough to resolve the 
dissipative scales), with the help of a 20 hot-wire probe 
anemometer, specifically designed for the particularities 
of the field experiment. 

Velocity gradients were computed without resort to the 
Taylor's frozen turbulence hypothesis. The Taylor-based 
Reynolds number of the flow, e stimated from the Tay- 
lor length A = y/uy{{diuiy) is Rx = 3.4 x 10^ (ui 
is the projection of the velocity fluctuations along the 
flow direction). We note that since the flow is somewhat 
anisotropic, the definition of a meaningful Taylor-based 
Reynolds number may be problematic. We will get back 
to this point in Sec. IV. 

The experimental velocity gradient pdfs are shown 
Fig.l. We find a good (within error bars) collapse of 
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FIG. 1: Experimental velocity gradient pdfs for atmospheric 
surface layer flow with R\ = 3.4 x 10^. Black line: sn; col- 
lored lines: Sij, with i ^ j. 
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FIG. 2: Numerical velocity gradient pdfs for homogeneous 
isotropic turbulence with R\ = 240. Black lines: diagonal 
components su', red lines: non-diagonal components Sij. 



standardized pdfs of velocity gradients Sij = djUi with 
i ^ j. Due to anisotropy effects in the surface layer, 
however, there is no collapse for the standardized pdfs 
of diagonal components, sa, and we have discarded the 
curves for S22 and S33 assuming, as a working hypothe- 
sis to be tested a posteriori, that isotropic results would 
correspond to the set {sn, Sy}, with i ^ j. 

The central aim of this work is to model the pdfs de- 
picted in Fig. 1 using direct numerical simulation (DNS) 
results for homogeneous and isotropic turbulence ob- 
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FIG. 3: Comparison between the numerical {Rx — 240; solid 
line) and the experimental {R\ — 3.4 x 10^; dots) pdfs of sn. 
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FIG. 4: Comparison between the numerical {Rx — 240; solid 
line) and the experimental {R\ — 3.4 x 10^; dots) pdfs of S23. 



tained at the considerably lower Taylor-based Reynolds 
number R\ = 240 (the numerical data corresponds to 
simulations discussed in Ref. [l3|)- The corresponding 
DNS velocity gradient pdfs are shown in Fig. 2. As it fol- 
lows from this figure, the pdfs collapse into two distinct 
groups, associated to the diagonal and non-diagonal com- 
ponents of the velocity gradient tensor Sij. Of course, we 
do not expect that the pdfs given in Fig. 2 yield a direct 
fitting to the ones of Fig. 1 - there is a clear discrepancy 
as shown in Figs. 3 and 4. 



IV. MONTE-CARLO PDF RECONSTRUCTION 

Our computational strategy is to consider the experi- 
mental = 3.4 X 10^) and the numerical {Rx = 240) 
flows discussed in Sec. Ill as the systems A and B, re- 
spectively, of Theorem 2. An important parameter here is 
the cascade distance Na — Nb of these flows. This quan- 
tity can be computed by measuring the flatness factors 
H4 of flows A and B and using them as input parameters 
in (I2.18p . From the pdfs of sn, we get 

H4{A) = 11.5 , 

i?4(S) = 6.6 . (4.1) 
Therefore, using (|2.18p . with = 1/3, we find 

iVA -iVs-^ logs (^^) =6.16. (4.2) 

Due to the discrete structure of the cascade in the multi- 
plicative models, we take 7V4 — Nb = 6 in the following 
considerations. 



1 




FIG. 5: The numerically reconstructed pdf of sn (black solid 
line) is compared to the experimental pdf (red circles). 

Using a random Poisson variable generator as the one 
given in Ref. [1^, it is straightforward to establish a 
stochastic process with random variable given by (I2.25P . 
On the other hand, in order to generate a stochastic pro- 
cess with random variable described by the numerical pdf 
of fill, we proceed in two steps: first, we define an accu- 
rate polynomial fitting to the logj^g profile; sec- 
ond, the polynomial analytical distribution just obtained 
is used in a Monte-Carlo accept-reject algorithm [25j . 
which produces random variables distributed according 
to Pb(sii). Analogous computations are performed for 
the numerical pdfs of S231 which are taken as a repre- 
sentative of the non-diagonal components of the velocity 
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FIG. 6: The numerically reconstructed pdf of S23 (black solid 
line) is compared to the experimental pdfs of Sij, with i ^ j 
(coUored symbols). 



gradient tensor. By multiplying the stochastic processes 
associated to the Poisson and the numerical pdfs we get 
standardized pdfs which would hopefully fit the experi- 
mental curves. We have taken a process with 2 x 10'' el- 
ements. In fact, an excellent agreement is attained from 
the Monte-Carlo reconstructed pdfs, as shown in Figs. 5 
and 6. A comparison between the modelled and the ex- 
perimental pdfs is also shown in Figs. 7 and 8 in linear 
scales, to be contrasted to Figs. 3 and 4. 
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FIG. 7: The numerically reconstructed pdf of sn (solid line) 
is compared to the experimental pdf (dots) in linear scales. 



FIG. 8: The numerically reconstructed pdf of S23 is compared 
to the experimental pdf (dots) in linear scales. 



It is important to emphasize that the remarkable fit- 
tings shown in Figs. 5-8, between the numerical and 
experimental pdfs for the set {sii,Sy}, are obtained 
from the mapping, determined by the single parameter 
Na — Nb, provided by the fluctuations given by (|2.25|) . 
This constitutes strong evidence for the existence of an 
underlying log-Poisson cascade process. We note, fur- 
thermore, that the agreement between modelled and ex- 
perimental pdfs would be not so good if the experimental 
pdfs of S22 or S33 were chosen in place of the one for sn. 
The present method, thus, has the heuristic potential to 
address issues of isotropy in boundary layer flows. 

A further application of our results is the definition of 
an effective Reynolds number Rx for the atmospheric sur- 
face turbulent flow, taking the more controlled Reynolds 
number of the DNS as a standard. We write, according 
to (ITTi)) . 

i?A = 240 xf— j -1.2x10^. (4.3) 

It was noted, in Ref. [l3|, that the rough estimate Rx — 
3.4 X 10^ displaces the point {Rx, H^) = (3.4 x 10^, 11.5) 
out of the empirical curve well modelled Hj^ ^ i?""* . How- 
ever, we find that if the alternative value (|4.3p is used 
instead of Rx = 3.4 x 10^, then the point {Rx,Hi) gets 
closer to the usual curve of flatness. 



V. CONCLUSIONS 

We have used the log-Poisson model of the turbulent 
cascade to get the pdfs of velocity gradient fluctuations 
of a high Reynolds turbulent atmospheric flow. The ex- 
cellent fittings are achieved by means of a Monte-Carlo 
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integration procedure and the use of standard pdfs ob- 
tained in a lower Reynolds number DNS. Our results indi- 
cate that non-gaussianity and anomalous scaling of scale 
dependent observables can be seen as different manifes- 
tations of intermittency that can be approached within 
a unified framework. Actually, this point of view has 
been formerly pursued along the multifractal description 
of intermittency (26j , with modest success in the quality 
of pdf fittings, nevertheless the fact that they are depen- 
dent on a large number of free parameters. 

As a natural application of our methodology, we have 
found a way to (i) select isotropy sectors of the veloc- 
ity gradient tensor in boundary layer flows and (ii) un- 
ambiguously define effective Taylor-based Reynolds num- 
bers in the presence of anisotropy. These results can be 
of considerable interest in the study of anisotropy effects 
in turbulent boundary layers. It is also likely that the 
same ideas can be extended to the case of free shear tur- 
bulence. 

An interesting question is how low can be the DNS 
Reynolds number, while still leading to good velocity 
gradient pdf fittings for higher Reynolds number fiows, 
along the lines discussed in Sec. IV. An investigation of 
this matter could throw some light on the problem of ex- 
tended self-similarity [l^]. Also, we wonder if correlation 
effects in the velocity gradient time series could be mod- 
elled in similar ways. A promising direction here would 
be to link the Fokker-Planck approach to turbulent time 
series [28] with the log-Poisson cascade model. 

It is clear that the multiplicative cascade picture is 
worth as a phenomenological construction if a consistent 



meaning can be given to concepts like the inertial range, 
local cascade, and the universality of velocity structure 
exponents. However, recent work on the scaling be- 
havior of velocity structure functions suggests that iner- 
tial and dissipative range fluctuations could be coupled 
in a bidirectional way. It has been found in 29] that 
the scaling exponents measured in the inertial range are 
changed if strong dissipative events are discarded in the 
averaging procedure, indicating a "flow of influence" from 
the small to the large scales. 

In order to address further related studies, we note that 
a possible solution to these puzzling observations, saving 
the essence of the multiplicative cascade phenomenology, 
would rely on the usual definition of the energy dissipa- 
tion rate as the local dissipation rate averaged over 
volumes with linear sizes of the order of £rn = L/a™". 
Since the energy dissipation rate is long-range correlated, 
it is likely that events which have strong local dissipation 
rates turn to be correlated to strong events in the above 
(inertial range averaged) sense. 
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